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Abstract 



We study the effect of adiabatically resonant conversion in galactic halos 
of neutrinos at the highest energies (~ 10^'^ - 10^^ eV), when the v source 
is in the center of a galaxy. Using the standard neutrino properties and the 
standard cosmological scenario for the hot dark part of matter, we find that 
interesting conversions may take place just for neutrino parameters relevant 
to the solar and atmospheric neutrino problem. The effect is due to the large 
enhancement in the v density in galactic halos and to the form of the effective 
matter potential both below and above the pole of the Z resonance. 
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It has recently been suggested [|T|-§[ that annihilation of extremely high energy (EHE) 
( > 10^^ eV) primary cosmic neutrinos on a background of clustered relic neutrinos may be 
able to account for the detected flux of cosmic rays of the highest energy ( > 10^° eV) [^Q. 
The Fly's Eye and AGASA ^ experiments confirmed the existence of proton-like cosmic 
rays with energies > 10^'' eV which, being well above Greisen-Zatsepin-Kuzmin cut- 
off, may create a rather serious problem in terms of a possible astrophysical explanation. 
In the above scenario this problem is elegantly circumvented by the fact that the 

highest energy cosmic rays may actually be generated very close to Earth, possibly in our 
hot-dark-matter galactic halo In order that this new mechanism for generating the 

highest energy cosmic rays has a chance to work, a relic electron-volt mass neutrino would 
be needed as the gravitational clustering of relic cosmic neutrinos is more efficient for heavier 
and slower species. If the Hubble and age parameter allow for a critical-density universe 
as suggested by inflation, then the mixed Cold + Hot Dark Matter model with m,y^ ~ m^^ ~ 
2.5 eV (giving Q,, ~ 0.2) was found to be consistent with all available observations. Indeed 
this model was shown to agree well with observations of cosmic background temperature 
anisotropy on large scales performed by the COBE satellite ||^ as well as with those on 
the number density of clusters [^]. Although in accordance with the theory of formation of 
cosmic structure, the scenario is beset by the problem 0] that one requires a very large energy 
generation rate of high energy neutrinos (which poses a difficult astrophysical problem) and 
the fact that Z, decay produces a large population of 7r° secondaries and hence we have 
more photons (from 7r° decay) than protons. 

The interaction of the EHE neutrino with the relic one had been studied first in [T^ 



and reconsidered afterwards in [0. It was shown ||TT| that neutrinos with the standard 
electro-weak interaction can travel through cosmic distances with a little hindrance, and 
only those neutrinos with energies close to the Z resonance value may be absorbed while 
crossing galactic halos, leading to a dip in the spectrum around Eres — 10^^ eV. This idea 
was then taken over in to show that secondaries (such as protons and photons) of the 
EHE-relic neutrino annihilation in our galactic halo may be responsible for the observed 
high energy cosmic ray energy spectrum. 

Here we examine whether the condition for adiabatically resonant conversion of these 
EHE neutrinos while traversing the halo around the source is met, if the oscillation param- 
eters are taken to be just those responsible for the solution to the solar and atmospheric 
neutrino anomaly. In the following we will be considering z/g — z^r(Ai), — {i^s being a 
sterile neutrino) as well as — Ug oscillations. 

Assuming that the halo is nearly round we take for the halo density a density profile like 

1 + [r/ay 

where a is the core radius of the halo. By considering the values for our local galactic halo 
it was found [Q that, notwithstanding the fact that light neutrinos constitute a smaller 



fraction of the density of the galactic halo |]T2| then they do of the cosmological density 



in a typical mixed-dark-matter scenario [13|, the neutrino number density inside the core 
might be a factor ~ 10^ — 10'' larger than the average cosmic v density. Hence, Ny ~ 
10^ — 10^ cm~^ was assumed in [Q, together with a relic neutrino mass = 10 eV. In 
the present work, however, we treat N^, as a free parameter, presenting first the conditions 
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under which conversion may occur, and then adding a discussion of whether or not these 
conditions may be met in 'reahstic neutrino halo' allowed by the Tremaine-Gunn phase-space 



constraint |]T4[. In addition, we assume that the halo consists predominantly of heavier ( 
and neutrinos with m,^^ ~ m^^ ~ 2.5 eV. It is just the fact that the background is not 
flavor symmetric that allows us to study not only active-sterile but active-active neutrino 
oscillations as well. 

Effects of a medium on neutrino propagation is determined by the difference of potentials, 
whose standard-model contribution in the context of Thermal Field Theory may readily 
be obtained from the relevant thermal self-energies of a neutrino: charged-current, neutral- 
current and tadpole. Of these, only neutral-current and tadpole contribution are of relevance 
here as the halo background consists of pure neutrinos (up to some of whatever dark-matter 
particles which appear as a major constituent of the halo and which are irrelevant for the 
present discussion). The contribution from the tadpole diagram, beside of being a constant 
independent of the external neutrino momentum and hence irrelevant (see below), is always 
proportional to the asymmetry for certain neutrino species. Since the cosmic neutrino 
background is likely to be CP-symmetric its contribution is completely negligible, leaving 
us with the neutral-current contribution only. 

Since the 'target' neutrinos are massive (m^^ ~ 2.5 eV) the Z boson may be considered 
"massless" always when s ~ 2Eymjj ^ M|. This corresponds to energies in the lab frame 
Ey > 10^^ eV. In the opposite limit, s -C M|, the Z boson should be considered "massive" 
and the usual contact approximation for the Z-propagator is adequate. Using the real- 
time version of Thermal Field Theory, we find by explicit calculation a contribution to the 
effective matter potentials in the "massless" limit as 

V; ^ 0.02 (Tj/E,) (^ |[/,,|2) ; (2) 

1=2,3 



V; ~ K ^ 0.02 (Tj/E,) . (3) 

Although the heaviest neutrinos are so massive that they are firmly non relativistic today 
they still should be described by the relativistic distribution function as their total number 
density is fixed at about 100 cm~^ in the uniform, non-clustered background and at about a 
factor of few larger in the galactic halo. Hence the temperature T^, from Eqs.(2) and (3) is 
always related to the neutrino number density via the relation, Ti, ~ (N^/O.ISY^^. Notice, in 
addition, that the appearance of a unitary mixing matrix in the vacuum U in the expression 
for Ve in (2) is due to the fact that internal lines of the thermal self-energy graph should 
necessarily be written in the (vacuum) mass eigenstate basis. If in a world with just three 
neutrinos all of them have nearly the same mass, of about ~ 1.5 eV, then the suppression 
due to {J2i l^^eiP) will be absent. 

Let us first consider the values of the parameters of the small mixing angle — Vs solution 



to the solar neutrino problem which lie in the region |15|: 



2.9 X 10"^ < /\mlJeV'^ < 7.7 x 10"^ ; 3.5 x 10"^ < sin^ 20 < 1.4 x 10'^ , (4) 

where Am^^ = — and 6 is the vacuum mixing angle of the z/g — z/^ mixing system. 
The MSW resonant condition for all neutrinos with > 10^^ eV now reads 
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0.04Tj(5:|f/e.P) ^AmL. (5) 

i 

Note that Eq.(5) is independent of neutrino energy. For the central value (5 x 10~^) in (4) 
we find that the MSW resonance for i/g -J^s oscillations occurs within the core radius of the 
halo if the neutrino number density inside the core of a galaxy is 3 x \Uei\^)~'^^'^ cm~^. 

The same occurs also for z/g — >■ oscillations as Vg > 0. 

The condition for unsuppressed oscillations is that the oscillation frequency must be real; 
otherwise the system is critically overdamped, the oscillations would be fully incoherent, 
and hence in fact there will be no oscillations. This translates into the requirement that the 
oscillation frequency should be much greater than the total neutrino interaction rate F, 

1 . . 



At the resonance this gives 



Am^ sin 26* 1 



For = 10^2 eV we have (Am^ sin 20)/2E^ ~ 2 x 10^^^ eV, whereas (1/2)F = (1/2) A^'^a ~ 
5 X (10"^2 - 10-3°) eV for ~ 10^ - 10^ cm'^ The latter estimate follows from p|, 
where the total neutrino absorption cross section was calculated to be cr{Ei, = 10^^ eV) ~ 
5 X 10"^^ cm^. Thus z/e — oscillations in the halo of a galaxy can be considered as 
unsuppressed for ~ 10^^ eV. However, level crossing alone is not sufficient for a 
complete conversion. The extent of conversion is determined by a degree of adiabacity; if 
there is only level crossing but no sufficient adiabacity at the resonance then there will be no 
complete conversion of electron neutrinos into sterile neutrinos. The condition for adiabatic 
transition at the resonance may be written in our case as 

3 Am2 sin^ 29 

where the subscript denotes that this quantity should be evaluated at the resonance. 
Using the distribution (1) one can easily show that for Ey ~ 10^^ eV adiabatically resonant 
conversion may only take place in the regions very close to the center of a galaxy, namely 
for 

ro < 2 X lO"-'^ a ~ 0(10 pc) . (9) 



For the last estimate in (9) we have taken a ~ 10 kpc [|T6|. Eq.(9) is due to the fact that the 



adiabacity parameter, as defined in (8), can be larger than unity for such EHEs only where 
the fiatness of the distribution (1) is appreciable. There is also another solution for ro from 
(8) but it is located well outside the core where the number density is smaller; besides, Eq.(l) 
is not a good fit to p(r) at r ^ a [0. However even the restriction (9) is consistent with the 



models ||T8| , p!9[l in which neutrino production takes place in Active Galactic Nuclei (AGN) 
jets. In these models protons are accelerated in the jets that stream highly coUimated out 
of their cores to distances of several parsec. In the jets the neutrinos are Lorentz boosted to 
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energies much higher than in AGN cores, what has important imphcations for their detection. 
Notice, for comparison, that K obtained here is always larger (for (Y,i \ Uei\'^Y^'^ ^ 10~^ and 
iVj, ~ lO'^ cm~^) than the effective matter potential in AGN cores whose magnitude was 
estimated to be 10~^^ — 10~^^ eV One should stress however here that as p{r) has a flat 
region for r < tq, the neutrinos are always produced at densities very close to the resonant 
one. Thus, the mixing angle at the production point 6mp — tt/4: rather than O^p ^ 7r/2, 
and consequently the oscillation probability P(z/e —>■ z/g) is centered around 1/2 rather than 
around sin^ 9. One can therefore speak of partial ( > 50%) conversion. 

Let us now turn to active-active z/g ^t(ii) oscillations. Notice that there is no resonant 
transition in the "massless" limit as Amg^^^-, > 0, but Vg — < [see Eqs.(2) and (3)]. 
The same holds for — ^ ^t{h) oscillations as Ve — V^(/i) < 0. The situation turns out to be 
completely different if we concentrate to the "massive" limit for the Z-boson propagator. 
In this case the effective matter potential is given as a sum of the finite-density and finite- 
temperature contribution 



where [for z^r(t'At)] 



and 



V{E,) = Vp + Vt{E,) , (10) 



Vp ~ V2GfN^L, (11) 



Vt{E,) ~ -2AGlTtE, . (12) 



In (11) A^^ represents the number density of photons and Ly the tau (muon)-like asymmetry. 
We are now in position to examine a possibility for resonant conversion as we have ~ 
~ -Vr > 0. It is to be noted that for = 10"^° and E^ ~ 10^° eV, \ Vt \, 
which represents a contribution coming from the low-energy tail of the Z resonance, is 
about 7 orders of magnitude larger than Vp for A^,^ ~ 10^ cm~'^ , and therefore from (10) 
V ~ 1/t < 0. 

Let us now concentrate to the hypothesis of two-neutrino vacuum oscillations which was 
found to provide a good quality description of the solar neutrino data for values of the two 



oscillation parameters belonging approximately to the region [^| 

5.0x10"" < Aml^^^-^/eV^ < 10"^° ; 0.65 < sin^ 2^ < 1.0. (13) 

The resonant energy for the central values in (13) (7.5 x 10~^^, sin^ 20 ~ 0.8) turns out 
to lie in the range 10^^ — 10^^ eV for the input neutrino number density 10^ — 10^ cm~^. 
The same occurs also for z/g u^i^i) oscillations as — Vf(^) > 0. Here we would like to 
stress the importance of Uf, — Uj. oscillations as the initial fluxes of the ultra-high energy 



neutrinos originating from AGNs are estimated to have a ratio [^ i^tI^h ^ 10~ (also 
Vf,/ Vp ~ 1/2 is expected from AGNs). Hence if an enhanced v-r/ ratio as compared to the 
non oscillation case is observed correlated to the direction of the source for ultra-high energy 
neutrinos, then this may represent a piece of evidence for neutrino oscillations taking place 



in the halo around the source. The model of Ref. [0 with maximum energy of 10 — 10 eV 



can therefore be very relevant to our scenario where level crossing occurs at a bit higher 
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energies. We find that resonant conversion with sufficient adiabacity does occur in the part 
of a dark halo with 



ro < 0(few pc) . (14) 

This is very similar to the estimate (9) for z/g — Us oscillations in the "massless" limit and 
hence still compatible to the proposed sites for neutrino production in the vicinity of AGNs. 
Let us finally consider the Vfj,—i's ||2^ oscillation hypothesis that provides an appealing ex- 



planation for the large up-down asymmetry observed by SuperKamiokande for atmospheric 
neutrino induced /i-events ||2^. The oscillations z/^ — Vg can solve the atmospheric problem 



if oscillation parameters lie in the approximate range 

3 X 10"^ < ^mlJeV'^ < 7 x 10"^ ; 0.82 < sin^ 2^ < 1.0 . (15) 

We note that it has been claimed |2^ recently that at present there is no constraint from 
cosmology to the z/^ — Vs solution with the aforementioned parameters. 

Since | 2E^V^ \«\ ^^fis 1 5 level crossing never occurs in a halo expect when sin^ 2^ is 
extremely close to the unity. Rather than considering such a coincidence we check out the 
adiabatic condition 

(sin^2^)-^^^^ » a-^. (16) 

While in the "massless" limit with E^, ~ 10^^ eV the condition (16) is well satisfied for the 
whole range < r < a , it is always extremely well satisfied in the "massive" limit. It 
is well known p6| that actually it does not matter whether a resonance has been crossed 
or not while a neutrino propagates adiabatically from one region to another as the similar 
formula for the propagation probabihty applies 

P{v^->y^,) = ^{1 + cos Oim cos 6 f ^) . (17) 

Eq.(17) describes a situation where z/^ produced at some density near the galactic center, 
with mixing angle 6im, propagates adiabatically through the halo to the region of diffuse 
neutrino background, with mixing angle 9fm- Since the mass squared difference Am^^ is 
always larger than the induced mass squared of the muon neutrino, the neutrino does not 
go through the resonance but is produced below it, and therefore we have 

e,„^ - Of^ ^9. (18) 

In this case the propagation probability (17) reduces just to the classical probability in the 
vacuum. With the parameters from (15) we get a constraint 

0.5 < P(z/^ ^ z/^) < 0.59 . (19) 

We have seen that all the above effects except that based on the z/^ — Ug oscillation 
hypothesis require at least N^, ~ 10'^ cm~^. There is, however, a strong bound based on 
Pauli's exclusion principle on the number density N'U^"'^ of neutrinos clustered in the halo 
around a galaxy. The neutrino density in the halo is approximately constrained because 



6 



of this requirement (Tremaine-Gunn constraint [T^) by N^""^ < {4:'K/3){pmax/h)^, where 



Pmax = ^vVesc IS the maximum momentum for a neutrino of mass rriy and v^sc is the escape 
velocity of bound neutrinos. If, in addition, we assume that massive neutrinos that fill up 
dark matter halos today result from gravitational instability of three cosmic seas of relic 
neutrinos, then the bound is improved by a factor of 2, giving 

^.™.<lxlO^(^)'(^)%n,-. (20) 

It is clear from Eq.(20) that to have sufficiently large Ny in halos where the neutrinos have 
typical escape velocities, larger values of uiu would be preferable. Conversely, if we stick 
to the Cold + Hot Dark Matter scenario of structure formation {my^ ~ rrtv^ ~ 2.5 eV), 
the increased neutrino density can be achievable only in halos where the neutrinos have 
large escape velocities. Notice, however, that larger values of my are still tolerable by the 
cosmological mass bound on light neutrinos (mj,^ + rriy^ + rriy^ < 37 eV), and thus for 
example rriy^ ~ rriy^ ~ 10 eV may result in Ny ^ 10^ cm~^. The effect of the higher 
neutrino masses would be to lower the resonant energy El'^^ (the energy which discriminate 
between the "massive" and "massless" limit) by a factor of few, thereby increasing the i^e — i^s 
conversion and depressing Ve — ^T{ti) conversion. 

In conclusion, we have considered a possibility for adiabatically resonant conversion of 
EHE neutrinos originating from the center of a galaxy, while traversing its dark matter halo. 
Owing to the large enhancement of the effective neutrino potentials at energies immediately 
above the pole of the Z resonance, we have found that z/g's can resonantly be converted 
into unobservable z/^'s, if their parameters are consistent with the MSW interpretation of 
the solar neutrino data. On the other hand, because the hot dark matter around galaxies 
is expected not to be flavor symmetric, z/g's can resonantly be converted into z^t-'s, if their 
parameters are consistent with the vacuum mixing solution of the solar neutrino anomaly. 
This kind of transition is especially important because there is a negligible number of z/t-'s in 
the initial fluxes of high energy neutrinos of astrophysical origin. Since the resonance density 
approaches the central neutrino production region, the above conversions turn out to be only 
partial. It is strikingly that in some proton blazar models a maximum neutrino energy is very 
close to the resonant energy. With regard to the Tremaine-Gunn constraint, these effects 
require either neutrino masses close to the cosmological mass bound or such halos where the 
neutrinos have large escape velocities, or both. Finally, if the parameters of active-sterile 
mixing system are consistent with the — Vg solution to the atmospheric neutrino anomaly, 
then neutrino propagation inside a dark matter halo is always adiabatic for Ey < 10^*^ eV 
and the propagation probability simply reduces to the classical probability in the vacuum. 

The author acknowledges the support of the Croatian Ministry of Science and Technology 
under the contract 1 - 03 - 068. 
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